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Abstract 

Purpose:lo assess the ability of a polarization transfer (PT) magnetic resonance spectroscopy (MRS) technique to improve 
the detection of the Individual phospholipid metabolites phosphochollne (PC), phosphoethanolamlne (PE), glyceropho- 
sphochollne (GPC), and glycerophosphoethanolamlne (GPE) in vivo In breast tumor xenografts. 

Materials and Methods: The adiabatic version of refocused insensitive nuclei enhanced by polarization transfer (BINEPT) 
MRS was tested at 9.4 TesIa In phantoms and animal models. BINEPT and pulse-acquire (PA) MRS was acquired 
consecutively from the same orthotopic MCF-7 (n = 10) and MDA-MB-231 (n = 10) breast tumor xenografts. After in vivo MRS 
measurements, animals were euthanized, tumors were extracted and high resolution (HR)-MRS was performed. Signal to 
noise ratios (SNRs) and metabolite ratios were compared for BINEPT and PA MRS, and were also measured and compared 
with that from HR-MRS. 

Results: BINEPT exclusively detected metabolites with ^H-^V coupling such as PC, PE, GPC, and GPE, thereby creating a 
significantly improved, flat baseline because overlapping resonances from Immobile and partly mobile phospholipids were 
removed without loss of sensitivity. GPE and GPC were more accurately detected by BINEPT in vivo, which enabled a reliable 
quantification of metabolite ratios such as PE/GPE and PC/GPC, which are important markers of tumor aggressiveness and 
treatment response. 

Co/7c/i/5/o/7; BINEPT Is advantageous over PA for detecting and quantifying the individual phospholipid metabolites PC, PE, 
GPC, and GPE in vivo at high magnetic field strength. As BINEPT can be used clinically, alterations in these phospholipid 
metabolites can be assessed in vivo for cancer diagnosis and treatment monitoring. 
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introduction 

Elevated concentrations of water-soluble phospholipid metab- 
olites such as the phosphomonoesters (PMEs) phosphochollne (PC) 
and phosphoethanolamlne (PE) and the phosphodiesters (PDEs) 
glycerophosphocholine (GPC) and glycerophosphoethanolamlne 
(GPE) are a metaboKc hallmark of cancer [1,2]. 'H and "P 
magnetic resonance spectroscopy (MRS) are able to detect this 
activated phospholipid metabolism in cancers in vivo, which can 
be used clinically for cancer diagnosis and treatment monitoring 
[2] . However, the total choline (tCho) signal detected by H MRS 
cannot be spectrally resolved into individual metabolites, mainly 
free choline (Cho), PC and GPC, due to the low spectral resolution 
in vivo at clinical field strengths of 1.5 and 3 Tesla (T), nor can it 



be resolved at high field strengths of 4 and 7 T [3,4,5]. Non- 
invasive detection of these individual phospholipid metabolites is 
of high interest as it has been shown that, for example, the PC/ 
GPC ratio can increase with increasing aggressiveness of breast 
cancer cells [6]. With '^'P MRS, individual PME and PDE signals 
can be detected, however, it suffers from relatively low signal to 
noise ratio (SNR) due to the intrinsically low sensitivity of the '^'P 
nucleus [7,8]. Even at high field strength, the detection of 
individual PE, PC, GPE and GPC is difficult, particularly in 
heterogeneous cancer tissues in which the homogeneity of die 
magnetic field is poor [9]. Another significant problem for 
quantitative in vivo detection with '^'P MRS is that these signals 
typically overlap with signals of other molecules with '^'P nuclei, 
such as sugar phosphates and immobile membrane phospholipids 
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(e.g. phosphatidylcholine, phosphatidylethanolamine), which cause 
a broad, uneven baseline, which in turn severely hampers the 
accuracy of metabolite measurement and quantification [10]. 
Accurate measurement and quantification of these phospholipid 
metabolites is particularly important in the clinic as consistent 
changes in phospholipid metabolite levels can aid in cancer 
diagnosis, prognosis, and treatment response monitoring. 

Another way to improve the SNR, baseline flatness, and 
detection accuracy of PMEs and PDEs is to apply polarization 
transfer (PT) methods [1 1] such as refocused insensitive nuclei 
enhanced by polarization transfer (RINEPT) [12,13]. PT 
techniques transfer the polarization of the excited 'H spins 
through J-coupling to the '^'P spins during the echo time period 
TE'H as shown in Fig. 1, which increases the SNR and removes 
all resonances without 'H-'^'P coupling, thereby flattening the 
baseline. Previous studies at lower field strengths have focused on 
the ability of 'H-^'P PT to increase the SNR. At 1.5T the SNR of 
PC obtained by using INEPT was comparable to that obtained by 
pulse-acquire (PA) acquisition [12]. However, for PDEs the SNR 
gain was 2 when using RINEPT [1 1]. At 3 T an SNR gain of 2.6, 
1.6, 2.2, and 2.4 for PE, PC, GPE, and GPC, respectively, was 
reported when comparing a selective-RINEPT to PA acquisition 
in human brain [10]. At 7 T the concept of PT was expanded by 
combining it with a direct '^'P detection sequence in one repetition 
time, which increases the SNR per unit of time as compared to a 
PT sequence alone [9]. 

As shown in the schematic overview of the adiabatic version of 
the refocused insensitive nuclei enhanced by polarization transfer 
(B INEPT) sequence in Fig. 1, the optimal duration for TE'H to 
achieve polarization transfer is relatively long due to the small J- 
coupling constants between '^'P and 'H spins, which causes 
substantial signal loss to molecules that have short T2 values for 
transverse relaxation. Consequently, apart from a higher SNR for 
PE, PC, GPE, and GPC, PT will simultaneously decrease the SNR 
of larger "P-containing molecules (with H-P coupling) that 
typically exhibit short T2 values and may overlap with these 
phospholipid metabolites [14]. This results in a flat baseline, which 
improves the detection of PMEs and PDEs. 

To test these considerations, we performed PT at the high field 
strength of 9.4 T on a preclinical scanner. We demonstrate that an 
optimized PT technique at 9.4T results in a flat and clean baseline 
without loss of sensitivity compared to Ernst-angle pulse-acquire 
acquisition, and in an improved in vivo detection and quantifi- 

BINEPT 



BIR4-1 




cation of PC, PE, GPE, and GPC in two breast cancer xenograft 
models. 

Materials and Methods 

Breast cancer xenograft models 

Highly aggressive MDA-MB-231 and weakly aggressive MCF-7 
breast cancer cells were purchased from the American Type 
Culture Collection and cultured as previously described [15]. 
2x10*^ MDA-MB-231 or MCF-7 cells in 50 nL of Hank's 
balanced salt solution (HBSS, Mediatech) were injected into the 
upper right thoracic mammary fat pad of anesthetized female 
athymic nude mice with an average weight of 26 g. Estrogen- 
dependent MCF-7 tumor growth was supported by a 0.18-mg 
17fi-estradiol 60-day release pellet (Innovative Research of 
America, Sarasota, FL) implanted subcutaneously into the back 
of mice 1 week before tumor inoculation. We studied 10 MCF-7 
tumors of 1.67±0.59 cc that developed within nine weeks and 10 
MDA-MB-231 tumors of 0.94±0.09 cc that developed within 
seven weeks. The maximum tumor size threshold for animals to be 
sacrificed was 2.5 cc. All surgical procedures and animal handling 
were approved by the Johns Hopkins University Institutional 
Animal Care and Use Committee, and conformed to the NIH 
Guide for the Care and Use of Laboratory Animals. 

Non-invasive in vivo IVlRS studies 

hi vivo '^'P MRS was performed on a 9.4T Bruker Horizontal 
Bore Small Animal Scanner (Bruker BioSpin Corp.). A 'H/'^'P 
double-tuned solenoid coil with an inner diameter of 12 mm was 
used (MRcoUs BV, Drunen, The Netherlands). The tumor was 
hanging into the coil while the animal lay on a platform above the 
coil with an opening for the tumor [16]. Mice were anesthetized 
by breathing isoflurane (2% in air) through a nose cone. Body 
temperature was maintained by using a blanket with circulating 
warm water. Breathing rate was monitored with a movement 
sensor attached to the abdomen of each mouse. 3D RARE images 
were acquired using an echo time (TE) of 7.2 ms, repetition time 
(TR) of 500 ms, RARE factor of 4, fUp angle of 90°, field of view 
(FOV) of 1 cnixl cnixl cm, 64 phase encode steps (64x64x64), 
and number of averages (NA) of 4. f and 2"** order Bq shimming 
was performed manually by minimizing the water line width of the 
entire tumor. Non-localized PA P MR spectra were acquired 
with adiabatic excitation (BIR4 45°, 200 |xs, 120 ppm bandwidth). 



BIR4-2/3 BIR4-4 BIR4 



31 



Figure 1. Schematic overview of the adiabatic version of the refocused insensitive nuclei enhanced by polarization transfer 
(BINEPT) sequence. Segmented BIR4 pulses (400 jis per segment, 35 ppm band widtfi) and a full BIR4 180" pulse (400 ^is, 35 ppm band width) 
were used instead of conventional 90 and 1 80 degree pulses. The echo time of 1 /4 J was set to 34 ms, which gave optimal signals for the acquisition 
of PE and PC. Abbreviations: BIR4, B1 insensitive refocusing; JpH, J-coupling constant; TE^H, effective echo time for proton refocusing; TE^V, effective 
echo time for phosphorus refocusing. 
doi:1 0.1 371/journal.pone.01 02256.g001 
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Figure 2. Raw RARE images and segmentation of corresponding slices. Necrotic tissue is depicted in red. Segmentation of (a) representative 
IVlCF-7 tumor and (b) representative iVlDA-IVlB-231 tumor, (c) Correlation betw/een line width of PCr and necrotic fraction for all MCF-7 and MDA-MB- 
231 tumors measured. 
doi:1 0.1 371 /journal.pone.01 02256.g002 



repetition time of 1 s, and 2000 averages, using a saturation slab to 
eliminate signals from adjacent muscles [17]. The combination of 
a short TR and 45° flip angle approximates PA acquisition with 
Ernst- angle excitation for PE and PC. Subsequently, an MR 
spectrum was acquired using an adiabatic version of the refocused 



insensitive nuclei enhanced polarization transfer technique 
BINEPT with a repetition time of 1 s and 2000 averages, which 
had an equal acquisition time of 33 min 20 s as PA. Segmented 
BIR4 pulses and a full BIR4 180° pulse were used as shown in 



PLOS ONE I www.plosone.org 



3 



July 2014 I Volume 9 | Issue 7 | e102256 



In Vivo 1H/31P Polarization Transfer iVlRS in Breast Cancer Models 

Fig. 1. Both TE'H and TE'^'P were set to 34 ms (Fig. 1), wliich 
was the optimal echo time for the detection of PC and PE [18]. 

Quantification of in vivo ^^P MRS data 

Lorentzian lines were fitted to the PA and BINEPT '"P MR 
spectra using JMRUI 4.0 software [19] and the AMARES 
algorithm [20]. For PA spectra, the phosphocreatine (PCr) 
resonance was set to 0 ppm, and line widths of PE, PC, GPE, 
and GPC were constrained to that of PCr during fitting (Fig. 2). In 
the BINEPT spectra, which only contained PE, PC, GPE, and 
GPC, line width was constrained to that of PE, and the chemical 
shifts of PE and GPC were estimated by JMRUI using a starting 
value of 6.8 ppm for PE and 2.9 ppm for GPC (Fig. 2). Frequency 
dilferences between the resonance positions of PC and PE, and 
GPC and GPE were fixed to 100 Hz in PA and BINEPT spectra. 
To be able to compare MR spectra of difiFerent mice with different 
coil loads and gain settings, metabolite levels were quantified as 
ratios to the noise, which was measured from the standard 
deviation of the last 200 points in the time domain signal. PE, PC, 
GPE, and GPC levels from the PA spectra were corrected for 
differences in Tl relaxation. P Tl values of PE, PC, GPE, and 
GPC were measured in vivo by progressive saturation series in 
four MDA-MB-231 tumors as shown in Table 1. PE, PC, GPE, 
and GPC levels from the BINEPT spectra were corrected for 
polarization transfer efficiency (about 50%) and for HTI, H T2, 
and ""P T2 relaxation with estimates based on literature values. 'H 
relaxation times of tCho were based on values obtained at 7 and 
9.4 T [21,22], and "P T2 relaxation times were based on values 
obtained at 3 and 7 T [10,23,24] (Table 1). 

Dual-phase extraction of tumors 

Directly after in vivo MRS measurements, mice were deeply 
euthanized as verified by toe pinch, and sacrificed by cervical 
dislocation. The entire tumor was immediately removed, freeze- 
clamped, carefully weighed for quantification, pulverized in hquid 
N2, and homogenized in 4 mL of ice-cold methanol with a tissue 
tearor homogenizer. Both lipid and water-soluble tumor extract 
fractions were obtained using a dual-phase extraction method 
based on methanol/chloroform/water (1:1:1; v/v/v) as previously 
described [25,26]. 

high-resolution (HR)-MRS studies 

The lyophilized, water-soluble extract fractions were dissolved 
in D2O containing 0.24x10 ' mol 3-(trimethylsilyl)propiomc- 
2,2,3,3-d4 acid (TSP; Sigma-Aldrich) and 6.0x10"** mol phenyl- 
phosphonic acid (PPA; Sigma-Aldrich) as internal reference 
concentration and chemical shift standards. ^'P HR-MRS was 
performed on a Bruker Avaiice 500 (11.7 T) spectrometer (Bruker 
BioSpin Corp.) using a 10-mm P probe. The MR spectra were 
acquired using the foUowuig acquisition parameters: 60° flip angle, 
10162 Hz sweep width, 8 K time domain size, repetition time of 3 
s and 2500 averages. The spectra were processed using the 
MestReC 4.9.9.6 software (Mestrelab Research). Lorentzian lines 
were fitted to the signals of PPA, PE, PC, GPE, and GPC. The 
areas under the curves were corrected for differences in Tl 
relaxation time and potential saturation effects. To this end, we 
measured the Tl relaxation times of PPA, PE, PC, and GPC with 
a progressive saturation series in a phantom solution in D2O. The 
Tl -corrected PE, PC, GPE and GPC levels were normalized to 
tumor weight. 
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Figure 3. Example of in vivo pulse-acquire (PA, top) and BINEPT (bottom) ^^P MR spectra of a representative IVICF-7 (left) and MDA- 
IVIB-231 (right) tumor. Lorentzian lines as fitted by JMRUI are shown below each MR spectrum. All phosphorylated metabolites are visible in the PA 
spectrum, whereas the BINEPT spectrum only contains signals from phospholipid metabolites with H-P-coupling such as PE, PC, GPE, and GPC. Note 
the broad, uneven baseline in the 0-5 ppm region of the PA spectra, where signals from mobile membrane phospholipids are resonating [35,45,46]. 
The signal of p-NTP is formed by p-NTP only. The signal labeled a-NTP is an overlapping signal from a-NTP, a-NDP, NAD, and DPDE. The signal labeled 
y-NTP is an overlapping signal from y-NTP and p-NDP. Typically, p-NTP is the smallest peak of the three NTP signals, however, here, y-NTP overlaps 
with a broad baseline signal that makes it appear smaller than p-NTP. 
doi:1 0.1 371/journal.pone.01 02256.g003 



MR data analysis 

To assess the quality of our in vivo MR data, we measured the 
line width of PCr in the in vivo PA data of each tumor. In 
addition, we calculated the necrotic fraction of each tumor by 
using the dark areas in the corresponding 3D RARE T 1 -weighted 
images of each tumor [27,28]. The tumor - air boundary was 
segmented in 3D using a threshold, set to the highest 10% in the 
histogram. Necrotic dark spots inside the tumor were segmented 
by manually adjusting a threshold to approximately the lowest 
10% in the histogram. Necrotic regions were quantified as necrotic 
fraction per tumor by counting the number of necrotic voxels in 
die dark spots inside the tumor divided by the total number of 
voxels inside the tumor. The correlation between PCr line width 
and necrotic fraction was measured for each tumor as shown in 
Figure 2. As the fitting of MR spectra becomes unreliable in cases 
where the line width is broad, the line width of PCr m the in vivo 
MR spectra was used as a criterion to exclude tumors with poor 
spectral resolution and large necrotic regions from further analysis. 
Data from tumors in which the line width of PCr in the in vivo PA 
MR spectrum was larger than 1 40 Hz were excluded from further 
analysis (vertical line in Fig. 2c), although when comparing the 
group results based on the data from all tumors, the conclusions do 
not change. This cutoff value in line width was a tradeoff between 
including as much data as possible and excluding data with poor 
quality. 

Statistical analysis 

A two-sided t-test assuming unequal variances between metab- 
olite levels or ratios was performed to test for statistically significant 
difierences between MDA-MB-231 and MCF-7 tumors, or PA 
and BINEPT methods. 



Results 

Our initial quality check of MR data demonstrated a strong 
positive correlation (Pearson correlation, r = 0.7053, P = 0.0011) 
between necrotic fraction and spectral line width as shown in 
Figure 2. This could be explained by the fact that in tumors with 
larger necrotic fractions more tissue-tissue transitions between 
viable and necrotic areas were present which degraded the Bo 
homogeneity, which in turn led to larger spectral line widths 
(Fig. 2). 

PA ^'P MRS, as shown in Figure 3, detected the PMEs PE and 
PC and the PDEs GPE and GPC, in addition to inorganic 
phosphate (Pi), PCr, and ot-, P-, and y-nucleoside triphosphates 
(NTP). Data from mice that did not pass the initial spectral quality 
check (line width of PCr<140 Hz) were removed from further 
quantification, which resulted in n = 6 for MDA-MB-231 and 
n=7 for MCF-7 tumors. BINEPT ""P MRS only detected PE, 
PC, GPE, and GPC. The PDF signals of GPE and GPC were 
better visible in the BINEPT spectra compared to the corre- 
sponding PA spectra as no overlapping signals from sugar 
phosphates, large membrane phospholipids, and other phosphor- 
ylated compounds were detected in the 0-10 ppm region, leading 
to a flat baseline throughout the spectrum (Fig. 2). Also, the PME 
in BINEPT spectra were not overlapping with other phospholipid 
signals. 

Figure 4 shows the quantification of the in vivo levels (average 
± standard error) of PE, PC, GPE, and GPC obtained with PA 
and BINEPT ^'P MRS. We did not detect any significant 
differences in PE, PC, GPE, or GPC levels in MDA-MB-231 
(n = 6) versus MCF-7 (n = 7) tumors. However, BINEPT '^'P MRS 
detected higher levels of GPE and GPC as compared to the 
detection by PA P MRS within the same tumor model. This was 
the case for both MBA-MB-23 1 as weU as MCF-7 tumors. The pH 
of these tumors as determined from the chemical shift difference 
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Figure 4. /n vivo phospholipid metabolite levels as measured by (a, c) PA and (b, d) BINEPT ^^P MRS. Average metabolite levels ± 
standard error are shown for (a, b) MCF-7 (n = 7) and (c, d) MDA-MB-231 (n = 6) breast tumor xenografts. Statistical results: PA vs BINEPT MCF-7: GPE 
p = 0.004 *, GPC p = 0.03 **; PA vs BINEPT MDA-MB-231: GPE p = 0.02 * GPC p = 0.06 **; MCF-7 vs MDA-MB-231 : PA, for all metabolites p>0.20; MCF- 
7 vs MDA-MB-231: BINEPT, for all metabolites p>0.21. 
doi:10.1371/journal.pone.0102256.g004 
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between PCr and Pi was 6.92+0.10 for MCF-7 tumors and 
7.00±0.11 for MDA-MB-231 tumors, which was not significantly 
different from each other. 



Figure 5 shows a representative P HR-MR spectrum of the 
aqueous phase of an MCF-7 breast tumor extract. The PME 
signals PE and PC as well as the PDE signals GPE and GPC were 
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Figure 5. ^'P High-Resolution iVlR spectrum of the same IVlCF-7 tumor as shown in Fig. 2. Thie full spectral range is displayed in the 
spectrum on tfie bottom. Signals of NTPs are visible in the 0 to -1 5 ppm region. An expanded region of 0 to 8 ppm is displayed on the top, as well as 
the region with the PPA reference signal at 13 ppm. The signals of PE, PC, GPE, and GPC are clearly separated. 
doi:1 0.1 371/journal.pone.01 02256.g005 



clearly separated in the P HR-MR spectra of tumor extracts 
(Fig. 5). Figure 6 shows the quantification of the PE, PC, GPE, and 
GPC concentration in MCF-7 and MDA-MB-231 tumors 
obtained from the '^'P HR-MR spectra. When normalized to 
tumor weight and the concentration reference PPA, no statistically 
significant differences in PE, PC, GPE, and GPC concentrations 
were detected between MBA-MB-231 (n = 6) and MCF-7 (n = 7) 
tumors (Fig. 6). 

Figure 7 shows a direct comparison between the in vivo PA, in 
vivo BINEPT, and tumor extract '^'P HR-MRS measurements, 
which were taken sequentially from the same tumors. The PE/ 
GPE ratio could not be assessed accurately in the PA spectra 
because GPE was barely detectable, restating in very small GPE 
values, leading to high ratios of PE/GPE and large standard 
dc\iations, ])articularly in MDA-MB-231 tumors. There was a 
difference between PC /GPC ratios obtained from in vivo 
BINEPT spectra and PA spectra. The value from BINEPT 
spectra more closely resembled the PC /GPC ratio that was 
obtained by HR-MRS of tumor extracts, for MCF-7 as well as 
MDA-MB-231 tumors. The PE/PC ratios obtained from in vivo 
PA and BINEPT measurements were similar. No significant 
differences were detected between any of the phospholipid 



metabolite ratios of MCF-7 versus MDA-MB-231 tumors. 
However there was a trend towards a higher PE/PC ratio in 
MCF-7 tumors as compared to MDA-MB-231 tumors. 

Discussion 

The polarization transfer technique BINEPT significantly 
improved the in vivo detection of the phospholipid metabolites 
PE, PC, GPE and GPC at high magnetic field strength in human 
breast tumor xenografts. The BINEPT sequence removed all 
signals originating from ""P compounds without P-H coupling 
because its detection is based on P-H polarization transfer. As a 
consequence, PCr, Pi, and NTP signals, as well as signals from 
sugar phosphates and large membrane phospholipids with broad 
line widths in the spectral region between 0 ppm and 10 ppm, 
were not detected by BINEPT, leading to a flat baseline, which 
significandy improved the quantification of PE, PC, GPE and 
GPC by enabling the use of simple Kne-fitting algorithms. 

An alternative way to obtain a flat baseline could be the use of 
an off resonance saturation pulse. The effect of such a pulse would 
be that magnetization transfer and chemical exchange would 
eliminate broad resonances from very large molecules [7]. 
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Figure 6. Quantification of PE, PC, GPE, and GPC levels from HR-MRS of (a) MCF-7 and (b) MDA-MB-231 breast tumor extracts. 

Average metabolite levels ± standard errors are shown and were expressed as concentration in ixmol per gram of tumor tissue. Statistical results: 
MCF-7 vs MDA-MB-231: PE p= 0.22 *, PC p = 0.77 GPE p = 0.31 *, GPC p = 0.16 *. 
doi:1 0.1 371/journal.pone.01 02256.g006 



However, the signal intensities of PE, PC, GPE, and GPC would 
be attenuated by such a pulse due to magnetization transfer, 
whereas BINEPT enhances these signal intensities. Proton 
decoupling applied during acquisition could improve the spectral 
resolution and SNR. However, at 9.4 Tesla the chemical shift 
anisotropy is the dominant relaxation mechanism, which wiU only 
result in a marginal increase of SNR, as the water line width is in 
the order of 80 Hz and the J-coupling in PMEs and PDEs is ~ 7 
Hz. In addition, the possible signal enhancement by nuclear 
overhauser enhancement (NOE) effect is different for each 
metabolite and spatially heterogeneous due to the use of a surface 
coU. This will introduce more variation into the data, which is not 
desired, and therefore proton decoupling and NOE were not 
applied in this study. Although the polarization transfer efficiency 
can also be different between compounds, this efficiency is easily 
calculated using spin operator calculations, and can therefore be 
taken into consideration during quantification. Only differences in 
T2 relaxation will attribute to an unknown variance in polariza- 
tion transfer efficiency. However, recent innovations in pulse 
sequence design may even recover this missing information in T2 
[23]. 

Considering absolute quantification, it is important to note that 
BINEPT, unlike PA, does not depend on the '^'P Tl relaxation, 
but rather on the T 1 relaxation of H and the T2 relaxation of P 
and 'H spins. We measured the "^'P Tl relaxation times at 9.4 T in 
vivo and estimated the ^'P T2 relaxation time to be 150 ms for PE 
and PC and 80 ms for GPE and GPC based on values reported in 
the brain at 3 T and the breast at 7 T [24,29] (Table 1). The 'H 
Tl relaxation time of the tCho resonance in the brain at 9.4 T is 



approximately 1400 ms [21] and its 'H T2 relaxation time at 7 T 
has been reported to be between 150 and 200 ms [22,30]. 
Therefore, we assumed the H Tl relaxation time of PE, PC, GPE 
and GPC in the tumor at 9.4 T to be 1400 ms, we estimated their 
T2 relaxation time to be 120 ms (shghdy lower than at 7 T). The 
fact that we used such estimates made it questionable to direcdy 
compare the quantified phospholipid metabolite levels from in 
vivo BINEPT, in vivo PA, and tumor extract HR ^'P MRS 
measurements. 

BINEPT ="P MRS detected higher levels of GPE and GPC as 
compared to the detection by PA ^'P MRS within the same tumor 
model. This underestimation of GPE and GPC levels by PA is 
most likely due to limitations of fitting with the JMRUI algorithm, 
which typically underestimates the fitting of signals with narrow 
hne width and low SNR, such as GPE and GPC, which are 
overlapping with other broad signals. The 'H and ^'P T2 
relaxation times used in the quantification are influencing the final 
metabolite levels measured from the BINEPT spectra displayed in 
Figure 4. The relaxation times of GPE and GPC are typically less 
accurately measured than the relaxation times of PE and PC, 
which have a much higher concentration in most tissues and hence 
appear with higher SNR. Therefore, the PE/GPE and PC/GPC 
ratios were measured more reliably by BINEPT than by PA ^'P 
MRS in VIVO. BINEPT '^'P MRS was best suited for the detection 
of the PE/GPE ratio since in contrast to the spectia obtained with 
PA, both PE and GPE appeared with sufficient SNR in the 
BINEPT spectintm. Accurate detection of PE/GPE and PC/GPC 
ratios is of great importance in studying cancer tissues, as these 
ratios are strongly associated with cancer aggressiveness [2,6]. 
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Figure 7. Comparison of pKiospKioiipid metabolite ratios between in vivo PA, in vivo BINEPT, and tumor extract HR ^^P MRS 
measurements of the same (a) MCF-7 and (b) MDA-MB-231 tumors. Values are given as average ± standard error. The PE/GPE ratios are not 
shown for in vivo PA measurements because the GPE resonance could not be reliably detected in several PA spectra due to low SNR. Statistical 
results: PA vs BINEPT: MCF-7 PC/GPC p = 0.02 *, all other ratios in MCF-7 and MDA-MB-231 p>0.12; PA vs tumor extract: MCF-7 PE/PC p = 0.04 *, PC/ 
GPC p = 0.05 **, all other ratios in MCF-7 and MDA-MB-231 p>0.21; BINEPT vs tumor extract: MCF-7 PE/PC p = 0.05 all other ratios in MCF-7 and 
MDA-MB-231 p>0.24. 
doi:1 0.1 371 /journal.pone.01 02256.g007 



Comparing ratios of metabolites in tumor extracts and in vivo 
measurements showed that the PC/GPC ratio obtained with 
BINEPT better resembled the tumor extract PC/GPC ratio than 
did the PC/GPC ratio obtained with PA. The quantification of 
BINEPT MR spectra using estimated values for relaxation times 
could be used in longitudinal studies on the effects of a given 
anticancer treatment on choline and ethanolamine metabolism. 

The PC/GPC ratios detected in this study are in good 
agreement with a recent study of Morse et al [31], in which 
MCF-7 tumors also displayed a slightly higher PC/GPC ratio 
compared to that of MDA-MB-231 tumors. These findings do not 
match the PC/GPC levels reported in the MCF-7 and MBA-MB- 
231 breast cancer cell lines grown in cell culture [6]. This is not 
surprising in light of the various differences that exist between cell 
lines in culture versus tumor xenografts. The PE/PC ratios in 
MCF-7 and MDA-MB-231 breast tumor xenografts in our study 
were comparabk- at around 1.5, which points towards a similar 
activation in choline and ethanolamine phospholipid metabolism 
in both tumor models, with a higher level of PE than PC. This is in 
good agreement with the finding that choline kinase alpha, which 
is overexpressed and activated in human breast cancers [32], has 
the ability to use ethanolamine as well as choline as substrates to 
produce PE and PC, respectively [33]. 

In early attempts with in vivo ^'P MRS of human cancers, a 
combined peak for PME and PDE was detected [34,35] . Some early 
studies demonstrated the potential of ^'P MRS for monitoring the 
effects of anticancer treatments [36,37] by evaluating the ratio of 



PME to PDE, and for assessing proliferation by evaluating the ratio 
of PME to Y-NTP [38]. However, at that time, the spectral 
resolution was insufficient for resolving the individual metabolites of 
PE, PC, GPE and GPC. With recent technological developments of 
improved human high-field MR scanners, coUs, and pulse 
sequences for "P MRS, resolving these individual metabolites has 
become feasible in the clinical setting [9,13]. Furthcrmon;, using 
human tissue biopsies and magic angle spinning "^'P HR-MRS, it 
has been possible to use phospholipid metabolites as markers to 
discriminate between tumor types and grades [39,40], evaluate 
response to treatment [36], predict response to treatment [41], and 
assess resection margins during breast cancer surgery [42] . 

As 'H MRS is more sensitive than ^'P MRS, it has been more 
widely used to study tCho in breast cancer. tCho has been used as 
a marker for malignancy and treatment response, however, the 
data from different groups show large variability and discrepancy 
between the findings [3]. Moreover, the H MRS detection of 
tCho in shrinking tumors often fails due to large signals from lipids 
[43] and poor magnetic field homogeneity (Bg), both of which are 
hampering accurate detection of tCho levels in breast lesions. This 
inaccuracy in tCho detection is particularly limiting for the 
assessment of tCho endpoints in longitudinal studies of response to 
treatment [3]. As ""P MRS does not pose the problem of large 
overlapping lipid signals and provides detailed information on 
anaboUc and catabohc metabolites in phosphohpid metabolism 
such as PE, PC, GPE, and GPC, it may be a better-suited 
technique to evaluate the role of chohne- as well as ethanolamine- 
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related metabolites in breast cancer. Even though the sensitivity of 
^'P MRS is relatively low, the PME levels in tumors can be very 
high, which allows for the possibility of assessing phosphoUpid 
metabolism even in small tumors. Furthermore, in the near future, 
better coils and improved pulse serjuences will become available 
for clinical use, which will increase the sensitivity of MRS. 
Finally, the abihty to detect the individual phospholipid metab- 
olites PE, PC, GPE, and GPG would significantly improve the 
diagnosis, treatment monitoring, and assessment of tumor 
response to therapy of breast tumors. 
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